Omicron BA.1 and BA.2 variants increase the interactions of SARS-CoV-2 spike glycoprotein with ACE2.

SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.

SARS-CoV-2 Delta Variant Decreases Nanobody Binding and ACE2 Blocking Effectivity.

The Delta variant spreads more rapidly than previous variants of SARS-CoV-2. This variant comprises several mutations on the receptor-binding domain (RBDDelta) of its spike glycoprotein, which binds to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. The RBD-PD interaction has been targeted by antibodies and nanobodies to prevent viral infection, but their effectiveness against the Delta variant remains unclear. Here, we investigated RBDDelta-PD interactions in the presence and absence of nanobodies H11-H4, H11-D4, and Ty1 by performing 21.8 µs of all-atom molecular dynamics simulations. Unbiased simulations revealed that Delta variant mutations strengthen RBD binding to ACE2 by increasing the hydrophobic interactions and salt bridge formation, but weaken interactions with H11-H4, H11-D4, and Ty1. Among these nanobodies H11-H4 and H11-D4 bind RBD without overlapping ACE2. They were unable to dislocate ACE2 from RBDDelta when bound side by side with ACE2 on RBD. Steered molecular dynamics simulations at comparable loading rates to high-speed atomic force microscopy (AFM) experiments estimated lower rupture forces of the nanobodies from RBDDelta compared to ACE2. Our results suggest that existing nanobodies are less effective to inhibit RBDDelta-PD interactions and a new generation of nanobodies is needed to neutralize the Delta variant.

Binding Mechanism of Neutralizing Nanobodies Targeting SARS-CoV-2 Spike Glycoprotein.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human cells upon binding of its spike (S) glycoproteins to ACE2 receptors. Several nanobodies neutralize SARS-CoV-2 infection by binding to the receptor-binding domain (RBD) of the S protein, but how their binding antagonizes S-ACE2 interactions is not well understood. Here, we identified interactions between the RBD and nanobodies H11-H4, H11-D4, and Ty1 by performing all-atom molecular dynamics simulations. H11-H4 and H11-D4 can bind to RBD without overlapping with ACE2. H11-H4, and to a lesser extent H11-D4, binding dislocates ACE2 from its binding site due to electrostatic repulsion. In comparison, Ty1 overlaps with ACE2 on RBD and has a similar binding strength to ACE2. Mutations in the Alpha variant of SARS-CoV-2 had a minor effect in RBD binding strengths of ACE2 and nanobodies, but reduced the ability of H11-H4 and H11-D4 to dislocate ACE2 from RBD. In comparison, the Beta variant weakened the RBD binding strengths of H11-H4 and H11-D4, which were less effective to dislocate ACE2 binding. Unexpectedly, mutations in Beta strengthened Ty1 binding to RBD, suggesting that this nanobody may be more effective to neutralize the Beta variant of SARS-CoV-2.

Critical Interactions Between the SARS-CoV-2 Spike Glycoprotein and the Human ACE2 Receptor.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects human cells by binding its spike (S) glycoproteins to angiotensin-converting enzyme 2 (ACE2) receptors and causes the coronavirus disease 2019 (COVID-19). Therapeutic approaches to prevent SARS-CoV-2 infection are mostly focused on blocking S-ACE2 binding, but critical residues that stabilize this interaction are not well understood. By performing all-atom molecular dynamics (MD) simulations, we identified an extended network of salt bridges, hydrophobic and electrostatic interactions, and hydrogen bonds between the receptor-binding domain (RBD) of the S protein and ACE2. Mutagenesis of these residues on the RBD was not sufficient to destabilize binding but reduced the average work to unbind the S protein from ACE2. In particular, the hydrophobic end of RBD serves as the main anchor site and is the last to unbind from ACE2 under force. We propose that blocking the hydrophobic surface of RBD via neutralizing antibodies could prove to be an effective strategy to inhibit S-ACE2 interactions.

Conformational transition of SARS-CoV-2 spike glycoprotein between its closed and open states.

In 2020, the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected millions of people worldwide and caused the coronavirus disease 2019 (COVID-19). Spike (S) glycoproteins on the viral membrane bind to ACE2 receptors on the host cell membrane and initiate fusion, and S protein is currently among the primary drug target to inhibit viral entry. The S protein can be in a receptor inaccessible (closed) or accessible (open) state based on down and up positions of its receptor-binding domain (RBD), respectively. However, conformational dynamics and the transition pathway between closed to open states remain unexplored. Here, we performed all-atom molecular dynamics (MD) simulations starting from closed and open states of the S protein trimer in the presence of explicit water and ions. MD simulations showed that RBD forms a higher number of interdomain interactions and exhibits lower mobility in its down position than its up position. MD simulations starting from intermediate conformations between the open and closed states indicated that RBD switches to the up position through a semi-open intermediate that potentially reduces the free energy barrier between the closed and open states. Free energy landscapes were constructed, and a minimum energy pathway connecting the closed and open states was proposed. Because RBD-ACE2 binding is compatible with the semi-open state, but not with the closed state of the S protein, we propose that the formation of the intermediate state is a prerequisite for the host cell recognition.